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EVOLUTION OF FLUORINE IN THE GALAXY WITH THE ^-PROCESS 
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ABSTRACT 

We calculate the evolution of fluorine in the solar neighborhood with the ^-process of core-collapse 
supernovae, the results of which are in good agreement with the observations of field stars. The v- 
process operating in supernovae causes the [F/0] ratio to plateau at [0/H] < — 1.2, followed by a rapid 

increase toward [O/H] ~ —0.5 from the contribution of Asymptotic Giant Branch stars. The plateau 
value of [F/O] depends on the neutrino luminosity released by core-collapse supernovae and may be 
constrained by using future observations of field stars at low metallicities. For globular clusters, the 
handful of [F/O] measurements suggest that the relative contribution from low- mass supernovae is 
smaller in these systems than in the field. 

Subject headings: Galaxy: abundances — Galaxy: evolution — stars: abundances — stars: AGB and 
post-AGB — supernovae: general 



1. INTRODUCTION 

Most of the energy from core-collapse supernovae is 
released as neutrinos and anti-neutrinos (> 10 53 erg). 
However, the interaction of the neutrinos with mat- 
ter and the effects on the nucleosynth esis have only 
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elements such as Fe and a elements, but it will increase 
those of some elements such as B, F, K, Sc, V, Mn, and 
Ti. 

Fluorine is an intriguing, though currently poorly stud- 
ied element. Most studies of F are from cool stars in 
which F measurements are only available from the HF 
molecule near 2.3 microns. Often, the F abundance 
comes from a single HF line. Asymptotic Giant Branch 
(AGB) stars a nd massive stars hav e both been suggested 
to produce F (porissen et al.lH~992T ) but the F production 
has only b e en confirmed for AGB stars (|Werner et all 
[20051 [2001 lAbia et al.l I20T0I [Milt IQtsuka et all 1201 U 
and references therein). Both in low- mass and massive 
stars, 19 F is produced by core and shell He-burning at 
T > 1.5 x 10 s K, but is destroyed by a-captures once 

the temperatures exceed ~ 2.5 x 10 8 K. In AGB stars, 
there is a primary component produced by the 18 0(n, 
7 ) 19 0(^-) 19 F reaction (IGallino et aljf2010h . which is in- 
cluded^ the n etwork used to compute the AGB yields 
()Karakasll20Tfjt) . AGB models with initial masses of 
~ 4 — 7Mq destroy F by proton captures that occur 
at the base of the convective envelope (hot bottom burn- 
ing). In AGB stars the production of F is highly mass 
dependent, where F production peaks at ~ 3M Q at solar 
metallicitv (ILugaro et alll2004D . 

iKobavashi et all ~ 2011bl ) showed that since the AGB 
mass range that produces F is 2 — 4M , this con- 
tribution is seen only at [Fe/H] >^ — 1.5 in Galac- 



tic ch emical evolution models (see also iTravagli o et all 
1999), and that the F production from AGB stars is 
not enough to explain the observations around [Fe/H] 
~ 0. In the other G alactic chemical evolution models, 
iTimmes et al.l (|1995l) showed that their F yields of core- 
collapse supernovae with the z/-process were not enough 
to meet the observations of stars. Massive stars evolving 
as Wolf-Rayet (WR) stars will experience very strong 
stellar winds, which may p revent the destruct ion of F 
(|Mevnet fe Arnouldl 12000) . iRenda et all (|2004l ) showed 
that [F/O] could be enhanced at [O/H] > - 0.2 with the 

WR yields of iMevnet fe Arnouldl (|2000D in addition to 
the AGB yields. However, the contribution by WR to F 
ma y be reduced by inc l uding rotation in the stellar mod- 
els (|Heger et all [2001 iPalacios et all l2005h . Therefore, 
we do not include the yields of WR stars in this paper. 

In this paper we show the effects of the ^-process in 
the Galactic chemical evolution of fluorine using latest 
yields of core-collapse supernovae and AGB stars. The F 
yields of SNe la are also included but are very small. In 
§2 we briefly describe our supernova models, where the 
details of the models will be described in Izutani et al. 
in preparation. In §3 we show the results of our chemical 
evolution models of the solar neighborhood including our 
nucleosynthesis yields with the ^-process. §4 denotes our 
conclusions and discussion. 

2. THE ^-PROCESS 

Although the cross sections of neutrino-nucleus reac- 
tions are small, a large flux of nutrinos is released when 
the core of a massive star collapses to form a neutron 
star. For this reason, the i/-process can have a signifi- 
cant effect on the nucleosynthesis of core-collapse super- 
novae. Neutrinos are emitted not only from a collapsing 
proto-neutron star but a lso from the innermost regio n 
just above a black hole ([Surman fe McLaughlin [20051) . 
We adopt the ^-process up to 80 Kr i n our nucleosyn- 
thesis calculations as in lYoshida et all (f2008h . both for 
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the cases of supernovae (SNe, the explosion energy of 
E = 10 51 erg) and hypernovae (HNe, E > 10 51 erg). 
The neutrino luminosity is assumed to be uniformly par- 
titioned among the neutrino flavors, and is assumed to 
decrease exponential ly in time with a timescale of 3 sec 
ijWooslev et al.lrT990l ) . The total neutrino energy is given 
by a free parameter and in this paper we present two 
cases with E v — 3 x 10 53 erg, which corresponds to the 
gravitational binding ener gy of a IAMq neutron star 
(jLattimer fe Prakashll2001l) . and 9 x 10 53 erg as the max- 
imum possible effect of the (/-process. The neutrino en- 
ergy spectra are assumed to be Fermi-Dirac distributions 
with zero chemical potentials. The temperatures of v^.t, 
v^r and v ei i> „ are set to be Ty = 6 M eV/fc and 4 MeV/fc, 
respectively (jRauscher et al.ll2002D . Note that the v- 
cross sections contain some uncertainties (|Heger et al.l 
f2005h . 

In a supernova, neutrinos interact with heavy elements 
through neutral-current reactions, and scatter off nuclei 
in or near their ground state, which lead to the excitation 
of particle unbound states that decay by neutron, proton, 
or a emission: 

(Z, A) + v -> (Z, A)* + v' -> (Z, A - 1) + n + v' (1) 

->{Z-l,A-l)+p + v\2) 
->(Z -2,A-A) + a + v{3) 

Charged-current reactions of v e or v e with heavy nuclei 
also play a role in producing new elements. These reac- 
tions correspond to the inverse processes of electron or 
positron captures. The new products in excited states 
emit 7-rays, neutron, proton, or a particles to decay to 
the ground state. The capture reactions of the protons 
and neutrons produced though these neutrino reactions 
also enhance the abundances of some elements. For most 
nuclei, neutral-current reactions are dominant because of 
the contribution from all flavors of neutrinos and higher 
temperature of i/^ T and £^ )T than that of v e and v e . 

We calculate the nucleosynthesis of core-collapse super- 
novae with progenitor masses of M — 15, 25, and 50M Q 
and initial metallicities of Z = 0, 0.004, and 0.02 for SNe 
and H Ne. The nuclear network includes 809 species up to 
121 Pd ()Izutani et al.|[2009t uzuFani fc UmedalboiOf ). The 
yields are calculated w ith the same assumptions as in 
IKobavashi et~"aT1 ((2006): for SNe, the mass-cut is set to 
meet the observed iron mass of 0.07M Q . For HNe, the 
explosion energy is set to be 10 x 10 51 and 40 x 10 51 
erg for 25 and 50M Q , respectively, and the parameters 
of mixing fallback models are determined to get [O/Fe] 
= 0.5. Although there may be diversity in the mixing- 
fallback process (as in the case of faint supernovae, e.g., 
IKobavashi et a l. 2011a), in this paper we focus on "typi- 
cal" supernovae that are dominant in the Galactic chem- 
ical evolution. 

In massive stars 19 F is mainly produced in a convective 
He shell as a secondary product through 15 N(a,7) 19 F, 
where the F yields are highly dependent on the metallic- 
ity. With the (/-process 19 F is produced in the O- and Ne- 
enriched region through 20 Ne((/, z/p) 19 F, and the F yield 
is increased by a factor of ~ 10 and 1000 for Z = 0.02 
and Z = 0, respectively. In the yields, the F/O ratio 
is smaller for more massive progenitors because of the 
larger mantle mass and larger O production, although 
the mass dependence of F/Fe is not so large. The F/O 



ratio does not strongly depend on the explosion energy, 
but F/Fe is smaller for HNe than SNe II because of the 
larger Fe production of HNe. 

3. GALACTIC CHEMICAL EVOLUTION 

We adopt the (/-process nucleosynthesis yields in the 
Galactic chemical evolution models. The nucleos ynthe - 
sis yields of AGB stars (1 - 7M Q ) from IKarakasI fl2010j ) 
are also included. We adopt the Kroupa initial mass 
function (I MF) and the same infall and star formation 
history as in IKobavashi et all ()2011bD . which reproduces 
the observed metallicity distribution function (MDF) in 
the solar neighborhood. 

Figure Q] shows the evolution of [F/O] against [O/H]. 
Without the AGB yields and the (/-process (short-dashed 
line), the predicted F abundance is too low to meet the 
observational data at all metallicities. With the AGB 
yields (long-dashed line), [F/O] shows a rapid increase 
from [O/H] > — 1.2 toward higher metallicities, which 

corresponds to the timescale of 2 — 4M Q stars in the so- 
lar neighborhood. At [O/H] - 0, [F/O] reaches -0.14, 
which is 0.26 dex larger than the case without the AGB 
yields. However, the present [F/O] ratio is still signifi- 
cantly lower than the observation s at [O/H] ~ 0. Note 
that compared to the yields from iKarakas fe Lattanziol 
(120071) . the F yields from AGB stars in IKarakasI (120101 ) 
were increas ed by applying the slower 19 F(a,p) 22 Ne re- 
action rate (jUgalde et all 12008! ). AGB stars may have 
polluted some Carbon-Enhanced Metal Poor (CEMP) 
stars with F at low metallicity via binary interactions 
(jLugaro et al.l 120081: iLucatelloetaLl 1201 ill , or through 
inhomogeneous enrichment. However, the overall con- 
tribution from AGB stars to the chemical evolution of 
the Galaxy is minimal at [Fe/H] > —1.5. 

The timescale of supernovae is much shorter than AGB 
stars, which means that the [F/O] ratio at low metal- 
licities can be strongly enhanced by the (/-process oc- 
curring in core-collapse supernovae. With the standard 
case of E v = 3 x 10 53 erg (solid line), the [F/O] ratio 
shows a plateau of [F/O] 0.4 at [O/H] < - 1.2, and 

reaches [F/O] ~ +0.19 at [O/H] > 0. This is consistent 
with the observational data of field stars at —0.5 < 
[O/H] < (ICunha et al.1 [20031: ICunha fe Smith! 120051: 
iCunha et al.l [2008) . If we adopt a larger neutrino lu- 
minosity of E v = 9 x 10 53 erg (dot-dashed line), [F/O] 
can be as large as ~ +0.37 at [O/H] ~ 0. 

In the bulge the star formation timescale is shorter and 
the average metallicity is higher than the solar neighbor- 
hood, but the [F/ O] ratio is not so differ ent at [O/H] 
~ (see Fig. 16 in IKobavashi et al.ll20lTbh . The obser- 
vations for the bulge stars (filled circles) might suggest 
that the IMF is also different, although the number of 
observations is too small to make a conclusion. 

At — 1 < [F/O] < — 0.5 the observational data are for 

stars in globular clusters (GCs), where the star formation 
and chemical enrichment histories are likely to be differ- 
ent to the solar neighborhood. These GC data seem to 
be more consistent with the models with the AGB yields 
only than with the (/-process. However, it is unlikely 
that the existence of the (/-process depends on the envi- 
ronment. With the (/-process the [F/O] ratio does not 
vary strongly with metallicity. Thus the differences ob- 
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Fig. 1. — Evolution of the [F/O] ratio against [O/H] for the solar 
neighbourhood with SNe II, HNe, and SNe la only (short-dashed 
lines), with AGB stars (long-dashed lines), with the v- process of 
SNe II and HNe (solid line and dot-dashed line for E v = 3 X 
10 53 and 9 X 10 53 erg, respectively). The dotted line is for the 
model for gl obular cluster s . The observational data sources are: 
open circles, ICunha et all J2003) . open squares, Cunha & Smith 
(2005) for the solar neighbourhood stars; filled circles. ICunha et alj 
(2008) for bulge stars. For the sta rs in globular clust ers, crosses, 
l Yong et al.l ||200ST) NGC 6712; plus. lSmith et al.l 12005T) . M4; stars, 
Cunha et al. (2003), <*> Cen; asterisks, Alves-Brito et al. (2011, in 
preparation), M22. 

served in [F/O] cannot be explained by variations in the 
metallicity of the progenitors. The neutrino luminosity 
may be small in the case of faint supernovae with a large 
black hole, which give high [a/Fe], but there is no signifi- 
cant difference seen in the [a/Fe] ratio between field halo 
stars and GC stars. One possible scenario is as follows: 
in GCs, the contribution from low-mass supernovae is 
smaller than in the field. Since the star formation occurs 
in a baryon dominated cloud with very high density, the 
initial star burst can be very intense. After the initial 
star burst, because of the small gravitational potential, 
outflow winds are generated immediately after the ex- 
plosion of massive supernovae, which may remove the 
contribution from low-mass supernovae. The small pro- 
duction of a elements from low mass supernovae means 
that the [F/O] ratio can reach values as large as ~ 0. In 
contrast, massive supernovae produce more a elements 
which results in [F/O] ratios of ~ —0.5, consistent with 
the observational data. 

The dotted line shows an example of such a GC model, 
where the timescale and duration of star formation is 
set to be 0.04 Gyr and 0.02 Gyr, respectively. This 
model gives the MDF peaked at [Fe/H] ~ -1.5. The 
[(Mg,Si,S,Ca)/Fe] is as large as in the solar neighborhood 
model, which is consistent with o bservations of GCs in 
our Galaxy (e.g- lPritzl et al.ll2005h . In a given GC, there 
is a spread in the observed O and F abundan ces, reflect- 
ing the so-called O-Na ant i-correlation (e.g.. iKraft et al.l 
U997tlCarretta et"aT1l2009[ ); there is a primordial popula- 
tion with high O and N along with low Na and Al, and a 



polluted population with low O and N along with high Na 
and Al. The polluted stars also include the products of 
H-burning at high temperature (~ 6.5 x 1 7 K), possibly 
from AGB stars or rotating massive stars (jGratton et al.1 
I2004D . Therefore, in Fig 1, we construct a model for 
GCs to fit the most O-rich stars, rather than the middle 
of the distribution. For the stars in M22 (asterisks) the 
pollution from AGB stars seems to be large. 

Note that the lack of low-mass supernovae is the oppo- 
site to the situation for dwarf spheroidal g alaxies (dSphs) 
which have low [a/Fe] and low [Mn/Fe] (fKobavash i et al.l 
2006). In dSphs, the dark matter component is large, 
the gas density is low, the star formation rate is low, 
and thus the contribution from massive supernovae is 
expected to be smaller than in the Milky Way halo. We 
do not include the pecul i ar sta rs with s-process contri- 
bution (jAbia et al.ll2010l 120111) and stars i n the Large 
Magellanic Could (e.g.. ICunha et all 120031 ) because in 
the first case F is produced by AGB stars, and in the 
second case the chemical evolution in the LMC may be 
quite different from the Milky Way. 

4. CONCLUSIONS AND DISCUSSION 

Both the ^-process of core-collapse supernovae and 
AGB stars play an important role in the production 
of fluorine. We succeed in reproducing the observed F 
abundances with our chemical evolution model that in- 
cludes the ^-process of E v = 3 x 10 53 erg. At low metal- 
licities ([O/H] < — 1.2) F production is dominated by 

supernovae, and thus future observations of field stars at 
low-metallicities are important for constraining the neu- 
trino luminosity released from a core-collapse supernova. 
If the neutrino luminosity is specified, the F abundance 
along with C could be a good clock in the study of galac- 
tic archaeology to distinguish the contribution from AGB 
stars and supernovae. The F observations of stars in GCs 
suggest that the star formation and chemical enrichment 
histories of GCs are different from those of field stars and 
that low-mass supernova played a smaller role in shaping 
the chemical evolution of these systems. 

The ^-process is also expected to be the producer of 
other elements such as K, Sc, and V. With E v — 9 x 10 53 
erg, [(K,Sc,V)/Fe] ratios are increased to be closer to the 
observational data, but such a large improvement is not 
seen with the standard value of the neutrino luminosity. 
There are several uncertainties that should be discussed; 
for K, the NLTE correction i n the observations is signif- 
icant ( Ko bavashi et al.l 120061) . The Sc yields could also 
be increased by the low-density mode ls that mimic 2D 
calculations (jUmeda &: Nomotol 120051 ) . There are also 
uncertainties in the reaction rates for V that may affect 
the nucleosynthesis calculations. 
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